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ABSTRACT 

We present the radial velocity confirmation of the extrasolar planet Kepler-447b, initially detected as a candidate by the Kepler mis¬ 
sion. In this work, we analyze its transit signal and the radial velocity data obtained with the Calar Alto Fiber-fed Echelle spectrograph 
(CAFE). By simultaneously modeling both datasets, we obtain the orbital and physical properties of the system. According to our 
results, Kepler-447b is a Jupiter-mass planet ( M p = 1.37^^ Mj up ), with an estimated radius of R p = 1.65^'j? Rj lip (uncertainties 
provided in this work are 3cr unless specified). This translates into a sub-Jupiter density. The planet revolves every ~ 7.8 days in a 
slightly eccentric orbit (e = 0.123^^3^) around a G8V star with detected activity in the Kepler light curve. Kepler-447b transits its 
host with a large impact parameter ( b = 1.076‘(° ( U?), being one of the few planetary grazing transits confirmed so far and the first 
in the Kepler large crop of exoplanets. We estimate that only around ~ 20% of the projected planet disk occults the stellar disk. The 
relatively large uncertainties in the planet radius are due to the large impact parameter and short duration of the transit. Planets with 
such an extremely large impact parameter can be used to detect and analyze interesting configurations such as additional perturbing 
bodies, stellar pulsations, rotation of a non-spherical planet, or polar spot-crossing events. All these scenarios would periodically 
modify the transit properties (depth, duration, and time of mid-transit), what could be detectable with sufficient accurate photometry. 
Short-cadence photometric data (at the 1 minute level) would help in the search for these exotic configurations in grazing planetary 
transits like that of Kepler-447b. 
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1. Introduction 

The unprecedented precision and long time span of the Kepler 
photometer have allowed the detection of exotic planetary sys¬ 
tems such as very_clqse-in hot-J upit ers around giant stars (e.g. 
Kepler-91b, iLillo-Box et al.l f2014cl lbh. super-Earth size plan¬ 
ets in the h abita ble zone of their parent stars (e.g., Kepler- 
22b. iBo rucki et ah 20121) . terrestrial worlds (e.g., Kepler- 
93b, Ballard et alj|20l4) . sub-Mercurv size planets (Kepler-37b, 
Barclay et al.ll^OOlT or very close-in packed planetary systems 


(e.g., Kepler-11, Lissauer et al.ll20lH) . This, combined with the 


more than 150000 stars that were continuously observed, favors 
the detection and characterization of planets in unusual configu¬ 
rations, difficult or impossible to detect with other techniques. 

In this context, the Kepler light-curve of Kepler-447 (KOI- 
1800, KIC 11017901; RA=19 /, 01 m 04'.46, DEC=48°33’36”) 
shows a V-shape d dip with a per iodicity of 7. 79430132 ± 
0.00000182 days (iBurke et alj|2014l) . According to lHuber et al.l 
(120141) . this object is a G8 Main-Sequence star slightly less mas¬ 
sive (Af* = 0.764 ^-q 45 M q ) and smaller (R* = 0.872^20 Rq) 
than the Sun. Usually, V-shaped eclipses are classified as false 
positives. They are mainly identified as eclipsing binaries with 
stars of similar sizes or grazing stellar eclipses. However, the 


* Based on observations collected at the German-Spanish Astronom¬ 
ical Center, Calar Alto, jointly operated by the Max- Planck-Institut fur 
Astronomie (Heidelberg) and the Instituto de Astroffsica de Andalucfa 
(IAA-CSIC, Granada). 


combination of transit and radial velocity measurements can re¬ 
veal the planetary nature of the transiting object. If a planetary 
mass is found for the transiting object, the V-shape is then ex¬ 
plained as a grazing planetary eclipse. To date, just one plane¬ 
tary grazing eclipse has been reported and confirmed to accom¬ 
plish the grazing criterion (h + RJR+ > 1) over 3cr (WASP-67b, 
iHellier etal .1120121 iMancini et al.l 20141) . Also, few other planets 
are known to transit their star in near-grazing orbits, but with 
the grazing criterion accomp lished below the 3cr level, as for in¬ 
stance WASP-34b (ISmallev et al.l|201 ill or HAT-P-27/WASP-40 
(iBekv et aPl201 lb I Anderson et al.l2oTlt) . So far, Kepler has dis¬ 
covered none. 

The discovery of such grazing planetary transits opens an 
interesting window. In this configuration, any (even small) grav¬ 
itational perturbation due to the presence of additional bodies 
in the system (like outer planets, exomoons, Trojans, etc.) would 
throw off the planet from its standard Keplerian orbit. Such effect 
would possibly induce a periodic variation in the impact parame¬ 
ter of the orbit that would be detectable in the transit data, so that 
they could become indirectly detectable (e.g., lKippindi2010l) . 
The observational imprints of other effects (such as planet rota¬ 
tion and oblateness or stellar pulsations) could be maximized in 
such interesting configurations. 

The present paper is part of an exoplanet survey car¬ 
ried out with instrumentation from the Calar Alto Observatory 
(CAHA). The follow-up of the Kepler candidates has already 
provided confirmation for two extrasolar planets by using the 
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CAFE instrument, namely Kepler-9 l b dLillo-Box et alJl2014hl) 
and Kepler-432b (ICiceri et alj 1201 5lP t has identified and an¬ 
alyzed some false positives and fast rotators dLillo-Box et al.l 
120151) : and has detected the secondary eclipse of WASP-lOb with 
OMEGA2000 (ICruz et al.ll2015l) . The present work, regarding 
Kepler-447, is organized as follows. In Sect. § [2] we summa¬ 
rize the ground- and space-based data used in this study. In Sect. 
§[3]we analyze the data by obtaining new spectrum-based stellar 
parameters (§ ED and modeling the radial velocity and transit 
signals in a simultaneous fit (§ 13.41 . In Sect. § 14. II we discuss 
the properties of this planetary system. In Sect. § 14.21 we discuss 
the possible (and exotic) configurations detectable in these ex¬ 
tremely grazing planetary transits and review an important sys¬ 
tematic effect causing the presence of outliers inside the region 
of short-duration transits. We summarize the conclusions of our 
analysis in Sect. §0 

2. Observations and data reduction 

2.1. Kepler photometry 

We have retrieved the Kepler photometry of this target from 
quarters Q1-Q17 (more than 1400 days, only with the long- 
cadence mode, 29.4 minutes). The publicly available data were 
downloaded from the Kepler MAST0 (Milkulski Archive for 
Space Telescopes). During Q4, the detector channels 5, 6, 7, 
and 8 in Module 3 of the charge-coupled device (CCD) failed. 
Consequently, any source falling in this detectors at any quar¬ 
ter could not be observed. Due to the rotation of the spacecraft, 
Kepler-447 felt in this broken down module on Q4, Q8, Q12, 
and Q16. Thus, no data were acquired for this KOI during those 
quarters. 

We used the simple aperture photometry (SAP) flux and 
its corresponding uncertainties provided by the Kepler team 
to compute the final light-curve (LC). We removed from this 
analysis any datapoint with a quality flag in the processed 
LC (SAP_QUALITY) equal to 128 (cosmic ray correction). In 
Sect. § 14. 3 1 we analyze the source and the effect of this data points 
on the LC. Artificial (well-known) large trends in the Kepler 
data and the observed modulations due to stellar activity were 
significant in this case. Consequently, we needed to detrend the 
fluxes to analyze the transit signal. We used a cubic spline func¬ 
tion to model the out-of-transit modulations, selecting the nodes 
by measuring the mean fluxes of one-day bins. This simple and 
quick approach provided a clean and flat light-curve in the out- 
of-transit region while keeping the transit unperturbed. 

2.2. High-spatial resolution imaging 

The host star candidate Kepler-447 was part of our high- 
resolution imaging survey of Kepler candidates (ILillo-Box et all 
l2012ll2014al) . We obtained lucky-imaging imagefjjfor more than 
170 KOIs with enough magnitude contrasts and spatial resolu¬ 
tion to statistically discard possible stellar configurations mim¬ 
icking a planetary transit. 

We defined the blended source confidence parameter (BSC) 
as the probability for a given object observed with high- 
resolution techniques and without any detected companions, 
to be actually isolated (ILillo-Box et al.ll2014al) . In the case of 
Kepler-447, the high-resolution images and their posterior analy- 

1 Also reported bv lOrtiz et all d2015l) 

2 https://archive.stsci.edu/kepler/data_search/search.php 

3 We used the AstraLux instrument at Calar Alto Observatory. 


sis revealed a BSC value of 99.9%. This implies a 0.1% of prob¬ 
ability of having a blended (undetected) source contaminating 
the Kepler light-curve. Thus, we will assume along this paper 
that Kepler-447 is an isolated target. 

2.3. High-resolution spectroscopy 

Thi s star wa s obs erved with the CAFE instrument 
dAceituno et al.l l2013t) installed at the 2.2m telescope in 
Calar Alto Observatory (Almerfa, Spain) as part of our follow¬ 
up of Kepler candidate^ This instrument has no movable pieces 
that could introduce artificial signals in the radial velocity. Also, 
it is located in an isolated chamber and its spectral coverage is 
fixed in the range ~ 4000-9500 A, with a spectral resolution 
in the range R = [59000,67000], depending on the order. The 
chamber is monitored in temperature, pressure, and humidity to 
check for possible relevant changes during the night. 

We obtained 21 high-resolution spectra of the target, sam¬ 
pling the different orbital phases of the planet candidate. The 
spectra were reduced by using the dedicated pipe line provided 
by the observatonQ. explained in lAceituno et alj (' 2013 b. We 
did not performed any dark correction and we used thorium- 
argon (ThAr) spectra obtained after every science spectrum to 
wavelength-calibrate the data. 

We cross-correlated the observed spectra against a weighted 
binar y mask to obta in the radial velocity (RV) of each epoch 
dBaranne et al.111996b . This mask was designed by using a solar 
spectrum obtained from the solar atlas provided by BASS 2OOC0 
(iDelbouillell 197 lh . From this spectrum, we selected a series of 
isolated, sharp, and high-contrast absorption lines. In total, our 
mask includes around 2100 lines in the CAFE spectral range 
(4000-9500 A). We used a velocity range of +30 krn/s around 
the expected radial velocity value for each case. The radial ve¬ 
locity of each spectrum is measured as the center of a gaussian 
fit to the resulting cross-correlation function (CCF), sum of all 
CCFs of the different orders. We corrected this velocity from the 
barycentric Earth radial velocity (BERV) obtained from the cor¬ 
responding Julian date at mid-exposure. In TableQ] we show the 
observing data and the final RV values for each epoch. 

3. Analysis 

3. 1. Stellar properties: rotation and physical parameters 

The SAP flux extracted by the Kepler pipeline shows a clear vari¬ 
ability at the level of ~ 4% (peak to peak, see Fig. Q}. These 
type of variations have been detected in many other sources with 
well-sampled photometric time series. They are related to the 
rotation of the star and the presence of stel lar spots on its sur¬ 
face. In particular, iMcOuillan et alJ ( 2013b ) analyzed the pho¬ 
tometric modulation caused by starspots in the Kepler sample 
of planet candidates. They calculated the rotational periods for 
several hundreds of KOIs by usi ng the autocorrelation fu nc- 
tion technique (ACF), described in IMcOuillan et all (l2013al) . In 
particular, they obtained a rotational period for Kepler-447 of 
P ro t = 6.459 ± 0.003 days. Added to this, we have computed 
the Lomb-Scargle periodogram of the light curve (see Fig. [2}. 
It shows a clear set of peaks around P rot , with the highest one 
at P peak = 6.4723 + 0.0003 days. Additional but less significant 

4 T his follow-up has alread y confirmed the plane tary nature of Kep ler- 
91 b dLillo-Box et all2014d lbl) and Kepler-432 b dCiceri et alj|2015l) . 

5 This pipeline was developed by Sebastian Sanchez. 

6 http://bass2000.obspm.fr/solar_spect.php 
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Fig. 1. Photometric time series of Kepler-447 as obtained by the Kepler telescope. The simple aperture photometry (SAP) is plotted in red. The 
quarters in which the source felt in the death detector channels of the CCD a re marked as shaded regions (see Sect. 12. It . The vertical dotted lines 
are plotted with a periodicity equal to the measured rotational period bv lMcOuillan et ah . J2013bl ). The orbital periods 55-59 are marked as blue 
filled circles (see Sect. 14.31 . 


peaks are found at larger periodicities, but their analysis is out 
of the scope of the current work. Both the orbital and the stellar 
rotation share similar periodicities. This could be due to some 
on-going synchronization between the stellar rotation and the 
orbital period due to similar processes studied in binary stars, 
such as tidal friction dZahnlfl966 ~) or hydrodynamical mecha¬ 
nisms (iTassoulll 1987b . 

At this point, it is worth to mention that the presence of 
stellar spots producing the high-amplitude photometric modu¬ 
lations implies the existence of stellar activity. This could con¬ 
taminate the transit signal with possible spot crossing events 
such as those detected in other hot Jupiters around active stars 
(e.g.. lSilva-Valio & LanzalfeOl ll : iDesert et al. 11201 ll) . These sig¬ 
nals can also be used to measu re sizes, temperatures, and posi¬ 
tions of the spots (lSilvall2003t) . as well as the spin-orbit align¬ 
ment between the planetary orbit and the stellar rotation axis 
(e.g., ISanchis-Oieda & Winnl20lTl) . However, the long-cadence 
data obtained for Kepler-447 and the short duration of its transit 
prevents this kind of studies in this system with the current data. 


Added to the stellar rotation, we used the CAFE spectra to 
obtain the physical properties of the host star. We combined 
all high-resolution spectra by shifting them according to their 
measured RV to obtain a moderate signal-to-noise (S/N) spec¬ 
trum. The final combination provides S/N = 65. This moder¬ 
ate S/N spectrum was used to obtain the stellar parameters of 
the host star. We derived the effective temperature, surface grav¬ 
ity, and metallicity by using s pectral syn th esis, using the Fe 
i spectral lines as described in IPene et al.! ( 2011 ). The results 
show that Kepler-447 is a G8V star with T e ff = 5493 + 62 K, 
logg = 4.40 ± 0.10, and [Fe/H] = +0.07 ± 0.05. These^values, 
agree with the photometrically obtained bv iHuber et al.l d2014l ) 
except for the metallicity, being sub-Solar in the latter work but 
with much larger uncertainties ([Fe/H]Hi 4 = -0.40/Jj ^). In Ta¬ 
bled we show both sets of stellar properties. 

By using these parameters, we can estimate the stella r mass 
and radius by applying the parametrization presented in Torresl 
d2010l) . We estimated the uncertainties in these parameters by 
running 10 5 Monte Carlo trials obtained from a gaussian distri¬ 
bution of the input values (i.e., T e g, log g, and [Fe/H]) and the 
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coefficients of the parametrization (see Table 1 in lTorresli2010l) . 
We then take the 68% limits (lcr) as the estimated uncertain¬ 
ties. According to this, we obtain M* = 1.00 ± 0.21 M 0 and 
Ri, = 1.03 + 0.16 R 0 . These values will be assumed along this 
work for further estimation of orbital and physical properties of 
the transiting companion. For completeness, we will also derive 
those properties by assuming the values provided bv lHuber et alJ 

( 12014 . 

3.2. Post-processing of the radial velocity and transit data 

CAFE data were collected during several runs in 2011, 2012, 
and 2013. Due to the long time span, we must account for sev¬ 
eral changes in the ThAr lamps along the whole set of obser¬ 
vations when comparing all RV epochs. In total, we can divide 
our data into three groups corresponding to three different ThAr 
lamps. We can correct the effect of the different lamps by apply¬ 
ing RV offsets between the samples. In total, we would need two 
offsets (v 0 ffi and v 0 ff 2 ) to set all three RV sets at the same level. 
These offsets will be included as free parameters in the RV fit, 
explained in the next section. 

By phase folding the detrended light-curve with the detected 
periodicity, we found few data points clearly over the transit 
signal. The source of these outliers (once the flagged data points 
with cosmic rays correction have been removed) is likely due to 
other instrumental effects. We ignored these flagged points when 
fitting the transit signal. Another remarkable feature to point out 
is that the transit is short (~ 1.13 hours) and V-shaped, rather 
than U-shaped as expected for planets. This suggests a large im¬ 
pact parameter, with the transit being grazing. 

3.3. Origin of the radial velocity signal 

We applied several line profile testfl to check for possible de¬ 
pendencies of the RV variation with stellar activity. We found 
no positive correlations between both sets of values in any of 
the performed tests. In particular, the bisector analysis (BIS, 
lOueloz et al.ll200ll ) provides a Pearson’s correlation coefficient 
of 0.007, thus showing no correlation of the RV with the line 
profile indicator. We show the BIS values against their corre¬ 
sponding RVs in Fig. [4] This may indicate that the RV is not 
correlated with possible line profile asymmetries caused by stel¬ 
lar inhomogeneities such as spots. 

We have also computed the periodogram of the radial veloc¬ 
ity data (see lower panel of Fig. S. It shows a relevant peak at 
the expected period from the transit analysis, with a false-alarm 
probability (FAP) of FAP < 0.01 %. This also indicates that the 
object transiting the host star is actually producing the RV vari¬ 
ations. 

3.4. Simultaneous fit of the radial velocity and transit signals 

Due to the aforementioned peculiar characteristics of this sys¬ 
tem, we decided to perform a simultaneous fit of the RV and 
transit data. This is because both effects share key parameters 
that could importantly affect the results of the other dataset (in 
particular the eccentricity and the argument of the periastron). 
The orbital period (P 01 b) and thejnid-transit time (To) were fixed 
to the values measured bv lBurke et al. ( 2014 ). 


7 By using the python code provided by P. Figueira (available at 
https://bitbucket.org/ pedrofigueira/line -profile-indicators), explained in 
iFigueira et ali ( 2013() and presented in ISantos et all (1 20141 ). 


In total, ten free parameters are needed to model both 
datasets: radial velocity semi-amplitude ( K ), eccentricity (e), ar¬ 
gument of periastron (to), systemic velocity of the system (Uys 13- 
two RV offsets (v' 0 ffi and v 0 ff 2 , see Sect. § ED>- semi-major axis 
to stellar radius (a/R*), planet-to-star radius ratio (R p /R*), in¬ 
clination (0, and a phase offset in the transit center to include 
possible unaccounted uncertainties in the determination of the 
transit epoch (c/> 0 ff). 

Among these free parameters , six are needed to model the 
transit signal by using the iMandel & Agoli (2002) formulation 
(i.e., Rp/R *, e, to, cilR+.j, and <poff)- According to the rec¬ 
ommendations from Muller et~akl (120131) for high impact pa¬ 
rameter transits, we fixed the limb-darkening coefficients to 
the theoretical values. We performed a trilinear interpolation 
of the determined stellar properties (T e ff, log g, and [Fe/H]) on 
the non-linear four term s limb-darkening values calculated by 
IClaret & Bloemenl (1201 ll) for the Kepler band. The resulting 
values used in the transit fit are a\ = 0.5895, «2 = -0.2477, 
«3 = 0.8538, and 04 - -0.4166. 

Due to the short duration of the transit (~ 1.13 hour) and the 
long cadence data obtained for this object (i.e., 29.4 minutes), 
there is a remarkable effect to take into account in the transit 
modeling. It is based in the fact that we try to fit long-cadence 
(timely binned) data with a theoretical transit model with in¬ 
finite resolution (i.e., unbinned). In practice, the long-cadence 
data smooths the transit signal, producing a broader transit. If the 
transit is grazing, this also translates intoa shallower depth. This 
effect was already highlighted by iKinnind (l2010h . who warned 
about large errors in the derived parameters if unaccounted. By 
definition, this binning effect is larger for short (with sampling 
step of the order of the transit duration ) and graz ing transits, such 
as the present case. As explained by iKinningl d2010l) . there are 
several ways to correct for it. The simplest (but most efficient) is 
to compute the light-curve model with a short cadence sampling 
(e.g., 1 minute) and then bin the model with the actual cadence 
of the observed data. We used this approach in our analysis. 



Fig. 2. Lomb-Scargle periodogram of the complete Kepler light curve 
(red line). A zoom into the shaded region around the rotational period 
of the star (—6.5) days is shown in the small panel. The vertical dashed 
line corresponds to the orbital pe riod, while the d otted ve rtical line is 
the measured rotational period bv lMcOuillan et all (2 013b ). 

Regarding the RV, we used the well-known model for the RV 
variations of a star surrounded by a single planet 

V(t) = V sys + K[cos (y(t) + to) + ecos w], (1) 


8 We note that this V sys is not the absolute systemic velocity of the 
system since it is not corrected from the possible instrumental offsets, 
although this should be of some tens of m/s. 
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Table 1. Measured radial velocities with CAFE. 


Julian Date 

S/N 

Phase 

RV 

Julian Date 

S/N 

Phase 

RV 

Julian Date 

S/N 

Phase 

RV 

(days)-2456000 



(km/s) 

(days)-2456000 



(km/s) 

(days)-2456000 



(km/s) 

76.593035 

26.0 

0.941 

1 421+ 0 - 027 
-0.027 

90.523633 

23.6 

0.729 


812.556916 

11.9 

0.365 

1 709 +0042 
i.z,wy _o.o42 

77.536041 

23.5 

0.062 

] t "3 c +0.032 

i . j jj_qq32 

91.376159 

18.8 

0.838 


817.599630 

7.0 

0.012 

1 301 +0056 

1.JV1_0 056 

80.543287 

25.7 

0.448 

1 361 +0031 

i.JOi_oo3i 

523.571588 

11.6 

0.288 

0 928 +0 043 

u.yz,o_0 043 

818.553484 

10.9 

0.134 

i i qq+0.034 

1 . 1 oO_o 034 

88.515082 

15.3 

0.471 

1 W +0 - 031 

17 -0.031 

597.302622 

18.2 

0.748 

l-136t«S 

818.585851 

11.4 

0.138 

1.1 OO_0 030 

88.549981 

17.4 

0.475 

1 389 +0033 

i.jo7_oo33 

801.517625 

9.5 

0.948 

1-308!°S 

823.582348 

13.8 

0.778 

1 462 +0025 

i.‘+uz,_0 025 

89.381498 

17.1 

0.582 

1 469+ 002 9 

1 .H-O7 _ 0029 

801.552053 

10.9 

0.953 

1 396 +0045 
i. jyo_oo45 

842.549205 

11.5 

0.213 

1 152 +0037 

1 . 1 JZ -0 037 

89.416506 

13.8 

0.587 

1 41 1+0.040 

1 11 —0.040 

812.524544 

11.9 

0.360 

1 303 +O033 

l.JUJ_0_033 

859.401764 

7.3 

0.375 

i i qo+0.039 
0.039 


where v(f) is the true anomaly of the planet. The true anomaly at 
each epoch was obtained by solving the Kepler equation accord¬ 
ing to the eccentricity and phase of the planet. 

We used our genetic algorithm GAbo$\ to explore the param¬ 
eter space and find the set of parameters that best reproduces the 
current data. GAbox allows for a clever exploration of this high¬ 
dimensional problem. We broadly restricted the parameter space 
according to our knowledge of the data. The adopted ranges for 
each of the free parameters are summarized in Tabled 

Due to the different number of data points on each data set 
(Nrv = 21 and Nlc = 994), we used the reduced chi-square as 
the minimization parameter in order to equally weight both ef¬ 
fects. Mathematically, this minimization parameter is expressed 
as Tred = -*id,Rv + ^ed.LC' The distribution of the 1500 con¬ 
vergence solutions found by GAboxin the parameter-parameter 
space is plotted in Fig. [5] In this figure we can see that all param¬ 
eters but the planet radius are very well constrained. However, as 
it is shown by the color-code in Fig.0 the least-square solution 
is located in the middle of the preferred region, suggesting that 
the adopted solution is a good compromise for all parameters 
OfredRV -1.2 and^ edLC = 1.8). We adopted this least-square 
set of parameters as our final solution. 

The uncertainties were obtained by running Monte-Carlo 
Markov Chain (MCMC) simulations with the Metropolis- 
Hasting algorithm. We ran 10 6 steps and discarded the 10% firsts 
to avoid dependence on the priors, which were set to the GAbox 
solution. We used 3cr levels as the adopted uncertainties. The fi¬ 
nal parameters are shown in Table[3] and the results of this fit are 
presented in Fig. [3] 

Given the obtained parameters in Table 0 we can derive the 
mass of the companion revolving around Kepler-447 by using its 
relation with the semi-amplitude of the RV signal 

f ,3 ~ n G ^sin 3 ' 

P(\-e 2 yr- (.Mp + M *) 2 ’ 1 J 

where the value inside the brackets is the so-called mass function 
(Fm). By assuming the stellar properties, we can obtain abso¬ 
lute values for the planetary mass and radius. Due to the relative 
d isagreeme nt between the stellar mass and radius obtained by 
iHuber et all (2 0141 ) and the present work (although still within 
the lcr uncertainties), we decided to get absolute parameters for 
both sets of stellar properties. However, we will refer to the fi¬ 
nal results as those obtained from our spectroscopic values (see 
Sect. § 13.11) . We thus obtain a mass of M p = 1.35^^ Mj up 
for the companion object, establishing its planetary nature. The 

9 GAbox is a fitting tool for any kind of scientific modeling 
and ha s already been used as spe ctral energy distribution fitter 
llRiviere-Marichalnr^tjilj lTQljiyTQjAl). transit and light-curve modu- 
lations fitting dLillo -Box et alj|2014cf). an d to extract the RV of high- 
resolution spectra iLillo-Box et alJ2014bl ). 


large uncertainty is consequence of the uncertainty in the stellar 
mass. Similarly, by using the R p /R * ratio, we obtain a planetary 
radius of R p = 1.65^ ^ Rj up . We should note that R p /R* is not 
very well constrained due to the short duration and grazing shape 
of the transit signal (see Fig. O. This is translated in a relatively 
large uncertainty for this parameter. 

4. Discussion 

4. 1. Derived properties of the planetary system 

We have shown that the companion transiting Kepler-447 has 
a planetary mass of M p - 1.37^ ^ Mj up and a large radius of 
R p = 1.65 1 d j 9 R Jup . According to these values, the resulting 
mean density of the planet would be p p = 0.30 pj up . This 
low density could indicate an inflated atmosphere for the planet. 
Although several inflation mechanisms could be taking place, 
the eccentricity found for this close-in giant plane t suggests that 
tidal heating could be the main mechanism (e.g., IJackson et al.l 
120081 ). Note that at periastron passage, the planet approaches its 
host star at a distance of r per /R* = 17.89^ ^. The non-negligible 
eccentricity found could be due to a third body in the system, 
preventing the circularization process and enhancing the tidal 
heating. 

The semi-major axis of the planetary orbit is a — 0.069+q^ 
AU as calculated by using the third Kepler law (assuming the 
mass of the host and the orbital period). We then obtain an im¬ 
pact parameter of b — 1.076^ ^, confirming that the transit is 
highly grazing. At mid-transit, only around 20% of the planet 
disk actually eclipses the star (i.e., A ec i/A p i = 0.20^j ^]). This 
explains the V-shape of the transit for this planetary mass com¬ 
panion. Given this impact parameter and the calculated planet- 
to-star radius ratio, we obtain a grazing criterion of b + Rp/R* = 
1.24^qJq (3cr uncertainties), clearly above the 3cr level. This con¬ 
firms that Kepler-447b is the second known planet with a grazing 
transit, and the first in the Kepler sample. The firstjplanet known 
to accomplish this criterion was WASP-67b (iHellier et al.ll2012l 
iMancini et al.l l2014h which was analyzed with short-cadence 
data from the ground. According to the Exoplanet ArchivcPl 
Kepler-447b would be the most grazing transiting exoplanet 
found to date, having the least fraction of its projected area cov¬ 
ering the stellar disk. We have illustrated this in Fig. [7] 

In Fig. [8] we show a scheme of the pole-on and face-on views 
of the system. As it is shown, the planet does not produce a sec¬ 
ondary eclipse due to the slightly inclined and the non-circular 
orbit. This is in agreement with the lack of a secondary eclipse 
for such a large and close-in planet. 

In Table [4] we present all derived physical and orbital prop¬ 
erties for this system by assuming both sets of stellar parameters 
from our own analysis of the high-resolution and moderate S/N 

10 http://exoplanet.eu 
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Fig. 3. Left panel: Phase-folded radial velocity data obtained during runs on 2012 (red symbols), 2013 (light blue), and 2014 (green). The black 
line shows the best fit model. The lower left panel shows the residuals of the model fit, having a rms of 31 m/s. Right panel: transit fitting of 
Kepler-447b. In the upper panel, we show the detected transit in the Kepler light-curve. We mark in color the outliers of the transit that have 
been removed from the fitting process (see Sect. 13.2k and that are mostly due to the misidentification of cosmic rays by the Kepler pipeline. 
The color-code is shown in the color bar and relates to the corresponding Barycentric Kepler Julian Date (BKJD = BJD-2454833 days). The final 
fitted model is represented by the solid black line. We have also included with a dashed blue line the original (not binned) model from which it is 
calculated (see Sect. 13.4b . The bottom panel shows the residuals of the fit, with a rms of 101 ppm. 


spectrum of the star (solution A) and from lHuber et al . ( 2014 1). 
solution B. 

4.2. Relevant configurations detectable in extremely grazing 
planetary transits 

The extremely large impact parameter that we have derived for 
the transit of the planet Kepler-447b means that interesting ef¬ 
fects non detectable in transiting planets with shorter impact pa¬ 
rameters can show up. In this section we summarize the most 
interesting effects that could be detectable in this kind of archi¬ 
tectures: 

1. Additional bodies- Due to the extremely high impact 
parameter of the orbit, small gravitational perturbations 
provoked by additional bodies in the system could induce 
the (projected) libration of the planet. This would induce 
periodic variations in the depth, duration, or transit time, 
depending on the relative position of both bodies. Among 
the different configurations with an additional third body in 
the system, the most common configuration could be the 
existence of an outer (non-transiting) planet in the system. 
In this case, we should also detect transit timing and depth 
variations, depending on the architecture of the system. 
However, short-cadence data is needed to analyze each 
transit individually and detect these variations. Another 
hypothesis can be the presence of a natural satellite orbiting 
around the planet (i.e., an ex omoon) that would induce an 
orbital wobble of the planet (iKinpin gl 120091) and thus the 
variation of the impact parameter on every orbit. In the 
most extreme case, the transit would not occur in some 
planetary revolutions due to this wobble of the orbital 


plane. Also, massive Trojan bodies can perturb the orbit 
of the plan et as it has been investigated, for instance, by 
lLaughlin & Chambersl (120021) . 

2. Stellar pulsation.- If the star is pulsating, the stellar radius 
will slightly increase and decrease, changing the R p /R * 
ratio and thus the fraction of the stellar disk occulted by 
the planet (or equivalently changing its impact parameter). 
These stellar pulsations should also be detectable in the 
out-of-transit time interval as periodic modulations of the 
light curve. 

3. Planet rotation and oblateness.- If the planet has a slightly 
ellipsoidal shape (i.e., non-negligible oblateness) and its 
self-rotation period is not coupled to the orbital period, the 
projected planet area could differ in the different transits. 
Thi s (som ehow exotic) idea has already been investigated by 
IZhu et all j2QL4j) fo r cand idates with short-cadence informa¬ 
tion. Also ICorreial (i2014l) studied the effects of non-spheric 
close-in planets on the transit signal. They suggested that the 
detection of the non-spherical shape could reveal the internal 
structure of the planet. This scenario would require several 
strong constraints on the planet rotation period. However, 
short-cadence (1 minute) precise photometric data is needed 
to measure the relevance of this effect. 

4. Starspots & spin-orbit orientation.- The occultation of 
starspots by the planet during its grazing transit increases the 
observed relative flux in the LC and could cause the pres¬ 
ence of brighter data points in the light curve. This scenario 
can be divided into two configurations: 
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- Spin-orbit alignment.- In the case of spin-orbit align¬ 
ment (where the orbit of the planet is perpendicular 
to the rotational axis of the star), the transit of the 
planet would only occult the higher latitudes of the star. 
This option would imply the existence of high latitude 
(polar) starspots remaining active during several months. 
Depending on the impact parameter of the transit we 
could measure the lowest latitudes covered by the 
planet. For instance, in the case of Kepler-447 the lowest 
region of the stellar disk covered by the planet would be 
^min = b - Rp/R * = 0.911, corresponding to a latitude 
of ~ 84°. Polar sta rspots have been found in main 
sequence stars (e.g.j_ Jeflfers^al]j2002^ and stars with 
rapi d rotation (ISch uessler & SolankiF 19921) . In particu¬ 
lar, ISanchis-Oieda et al.l ( 2013 ) detected spot-crossing 
events in the transits of the hot Jupiter Kepler-63b. 
The almost perpendicular orientation of this orbit with 
respect to the spin axis of the star indicated the presence 
of a long-lived starspot at one of the star rotation poles. 
The existence and persistence of these polar spots in 
Solar analogs has been pr evious ly theorized based on 
simulations (e.g., iBrown et al.l l20ld and references 
there in) and have been detected in few ca ses using 
different techniques (IStrassmeier & Riceii 19981) . 


- Spin-orbit misalignment.- The misalignment of the 
stellar rotation axis and the orbital angular momen¬ 
tum ha s be en obser ved in planetary systems (e.g., 
ISanchis-Oieda & Winn 20111) and particularly in hot - 
Jupiters (e.g.. iHebrard et al.ll2008t iTriaud et al.ll2010l) . 
A significant portion of these detections (~ 4 0%) re¬ 
sulted in obliquities larger than i/j ~ 22.5° (lAddison et al.l 
2014 ). In such case, the grazing planetary transit can 
shade lower latitudes of the stellar disk, where stellar 
spots usually exist. The imprints of this effect on the ra¬ 


dial velocity would be the presence of as ymm e tries in 
the Rossiter-McLaughlin effect (see, e.g., lOueloz et al.l 
l2000f) . If spot-crossing events are found in a planetary 
grazing transit this possibility can be tested by mea¬ 
suring the obliquity of the system with more accurate 
and higher-cadence RV measurements during the transit. 
However, we note that, as stated by IWinn et akl (12010l) . 
the great majority of planets around stars cooler than 
T etf < 6250 K are well aligned. 


Due to the short duration of the transit (7j ur = 1.135 + 0.016 
hours) and the long cadence observations for Kepler-447 (one 
data point every 29.4 minutes), it is not possible to perform 
an individual analysis of the different transits and thus test any 
of these interesting configurations. Short-cadence photometry 
would, instead, allow an individual fitting of the single transits. 


4.3. The effect of cosmic rays in short-duration transits 

The Kepler pipeline processes the data and tries to correct for the 
presence of possible cosmic rays inside the aperture of a target. 
This correction is sometimes affected by erroneous identification 
of cosmic rays in the transit time interval of long-cadence data 
when the duration of the transit is short as compared to the ca¬ 
dence. When this situation happens, the pipeline tries to correct 
the photometry and erroneously assigns larger fluxes for the in¬ 
transit data in the mentioned configurations. The result of this in 
the phase-folded light curve is a transit somehow filled by out¬ 
lier data points symmetrically distributed around the mid-transit 


time. These outliers in the transit region could be erroneously in¬ 
terpreted as hints for the detection of some of the configurations 
described in Sect. § 14.21 The Kepler team provides information 
about which data points have been corrected for the possible ex¬ 
istence of a cosmic ray by flagging them in the SAP_QUALITY 
field of the delivered data files with the flag 1210. 

Due to the very specific situations in which this effect can 
play an important role and contaminate the transit region of a 
planetary transit, there are few cases where this instrumental ef¬ 
fect has been reported. In particular, in lHerrero et al.l (i2013l) . the 
authors misidentify the erroneously corrected data points in the 
eclipse interval with a spot-crossing event in the system LHS 
6343. In a recent work bv lMontet et al.l (120141) . the authors iden¬ 
tified the data points corresponding to this spot-crossing event 
in LHS 6343 as flagged data points with the SAP_QUALITY 
equal to 128, and thus being cosmic rays misidentifications by 
the Kepler pipeline. 

In the case of Kepler-447, the version of the pipeline used in 
the first twelve quarters produced this kind of misidentifications, 
while no erroneous detection was found in the remaining quar¬ 
ters. The ~12% of the data points inside the transit time interval 
along the whole Kepler mission where tentatively identified (and 
thus corrected by the pipeline) as cosmic rays (see Fig. [0. The 
result, if the quality flags are not taken into account, is a tran¬ 
sit somehow filled by outlier data points (see Fig. [3] right panel) 
that could be erroneously interpreted as true signals caused by 
the configurations analyzed in Sect. § 14.21 Besides, owing to the 
extremely large impact parameter of the transit of Kepler-447b, 
thus being a planetary grazing transit, these flagged data points 
could by chance be mimicking the signal induced by some of 
those configurations. It is thus worthy to warn about this effect 
and to highlight the need of taking into account the quality flags 
provided by the Kepler team in short-duration transit events. 


5. Summary and conclusions 

We have confirmed the planetary nature of Kepler-447b by an¬ 
alyzing its radial velocity imprint on its G8V host star and the 
periodic dips observed in the Kepler light curve. Our simul¬ 
taneous modeling of both signals determined that Kepler-447b 
is a giant planet with a mass of M p - 1.35^^ and a radius 
of R p — 1.65^2jg, orbiting in a slightly inclined orbit. This, 
combined to its closeness to the host star ( a/R+ ~ 20.4), im¬ 
plies an extremely grazing transit, with a grazing criterion of 
b + R P /R * > 1, above the 3<x level. This is thus the most graz¬ 
ing transit known so far and the second confirmed plane t to ac - 
complish this criterion (after WASP-67 b, iMancini et alj|201~4h . 
Thus, we have found and characterized the first planetary graz¬ 
ing transit in the very large crop produced by Kepler. 

The extremely large impact parameter found for this planet 
can help in the search for additional effects playing a role in 
the system. We have discussed four interesting configurations 
that could be detectable by obtaining high-accurate and short- 
cadence photometry of grazing planetary transits. These include 
additional perturbing bodies (such as outer planets, exomoons, 
Trojan bodies, etc.), stellar pulsations changing the stellar radius 
and thus the impact parameter, planet rotation and oblateness, 
or spot crossing events. The non-negligible eccentricity found 
for this close-in giant planet suggests the presence of additional 
companion/s in the system. Thus, future works obtaining short 

11 See description of all SAP_QUALITY flags in 

the latest version of the Kepler archive manual in 

https ://archive. stsci ,edu/kepler/manuals/archive_manual .pdf 
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cadence photometry with the typical Kepler accuracy and during 
a relatively long time span would allow individual transit analy¬ 
sis and look for transit depth and duration variations as well as 
possible spot crossing events. 
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Table 2. Stellar parameters derived by different techniques. 


Parameter 

Huber2014 

This work ($13.11 

Feff (K) 

5555+J” 

5493 : 

t 62 

log(g) (cgs) 

4 440 + 0.119 

i +. t + i +u_o 308 

4.40 ± 

0.10 

[Fe/H] 

-0 40 +0 - 36 
-0.26 

0.07 ± 

0.05 

M* (M 0 ) 

0 - ? 64 

1.00 ± 

0.21 

R* (Re) 

0.872+O”« 

1.03 ± 

0.16 


Table 3. Parameters of the joint fitting of the radial velocity and transit 
data. 


Parameter 

Range 

Value 11 

Units 

K 

[50,250] 

141+ 42 

A ^ 1 -42 

m/s 

e 

[0.0,0.5] 

0 123 +0037 

-0.036 


CO 

[0,360] 

98 3 +u 

370 -11.0 

deg. 

Vsys 

[1.0,2.0] 

1 QQ+0.11 

km/s 

a/R * 

[10,40] 

20.41+jj-jg 


Voffl 

[-1000, 1000] 

-363+2^ 

m/s 

Vofi2 

[-1000, 1000] 

-84^ 

m/s 

R P /R * 

[0.01,0.50] 

O 

Os 

1 + 
o o 
© © 

00 >o 


i 

[80,90] 

86.55+Q32 

deg. 

0off 

[10,20] 

i 99+0.37 

A * -0.11 

min. 


Notes. (l,) Uncertainties are 3 cr. 

Table 4. Derived physical and orbital parameters from the fitted solu- 


tion of the joint analysis. 

Parameter 

Derived parameters 

Units 

Assumptions^- 1 


Sol. A ,a> 

Sol. B ,b > 



R* 

1 05 +019 

A U -0.19 

o 915+0.143 

Re 

M*, logg 

a 

0 0769 +0 0062 
u - 00079 

0 0703 + ° 0036 

u,u/uj -0. 00 23 

AU 

Porb, M* 

b 

1 076 +0112 
a.u /u_o 086 

1 076 +0112 
1.VJ/U_0 086 


a/R+, i, e, a) 

Tim 

i i -ic+0.016 

A -0.016 

i i q^+0.016 
A - A ' :,J -0.016 

hr 

TR model 

M„ 

1 ^7+0.48 

A * -0.46 

1 14+0.38 

1,1 -0.34 

^Jup 

K , M*, P mb 

Rp 

1 65 +0 59 

A-U -0.56 

1 40 +0 - 89 

A -^ -0.42 

^Jup 

Rp/ 

<p P > 

0 30 +0 71 

U J -0.24 

0 42 +0 - 74 

-0.37 

Ptup 

M p , Rp 

-^ecl/^pj 

20 +3 ° 

-23 

20 +30 
^ -23 

% 

b, Rp/R* 


Notes. All uncertainties cover the 99.7% of the probability distribution 
(i.e., 3cr). (a) Stellar radius and mass obtained from our spectroscopic 
analy sis (se e Sect . § 13.1b . <6) Stellar radius and mass obtained from 
iHuber et all d2014t> . (c) Assumed parameters in the estimation of the de¬ 
rived properties, t Fraction of the planet’s projected area eclipsing the 
stellar disk at mid-transit (in %). 
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Fig. 4. Upper panel: bisector analysis as a function of the measured radial velocity obtained during runs on 2012 (red symbols), 2013 (light blue), 
and 2014 (green). Lower panel: periodogram of the radial velocity data showing the significant peak at the corresponding transit period (vertical 
dashed line). 
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Fig. 5. Parameter-parameter representation of the 1500 convergence solutions of the transit fitting using GAbox. In color-code we represent the 
combined reduced chi-square statistic of the RV and transit fitting models (dark colors represent poorer chi-square models). The diagonal panels 
show the histograms of the distribution of the individual parameters within all solutions. The white circle represents the least-square model. 
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Fig. 6. Flux difference between the observed measurement and the transit model for data points inside the transit region. Outliers mostly due to 
the misidentification of cosmic rays by the Kepler pipeline are marked in the same color code as in the right panel of Fig. [3] The Kepler quarters 
of the mission with few or no data for this KOI due to the failure of CCD Module 3 are highlighted as gray shaded regions (see Sect. §ED- 
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Fig. 7. Left panel: impact parameter versus planetary density of known exoplanets with available information. Right panel: fraction of the 
projected planet transiting its host star as a function of the impact parameter. In both panels, Kepler-447b has been highlighted in red color and big 
circle symbol, including uncertainties. For reference, WASP-67b is also highlighted in light blue color with a square symbol. For clarity reasons, 
we do not include uncertainties for the other planets. 



Fig. 8. Orbital scheme the solution found for the orbit of the planet Kepler-447b (red filled circle). Left panel: Pole-on view of the orbit (separation 
from the star in units of the stellar radius). We show 500 orbits by bootstrapping the parameters inside their uncertainty limits. Right panel. Face 
on view of the orbit as projected in the sky (separation from the star in units of the stellar radius). We have marked the path of the planet between 
the dotted lines. The gray solid lines represent 500 orbits by bootstrapping the impact parameter. We can see that the transit is always grazing 
independently of the parameters taken inside the uncertainty limits. 
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